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In the present work, the masses of the decuplet baryons at finite temperature are 
investigated using thermal QCD sum rules. Making use of the quark propagator at 
finite temperature, we calculate the spectral functions to T 8 order, and find that 
there are no contributions to the spectral functions at T 8 order and the temperature 
corrections mainly come from that containing T 4 ones. The calculations show very 
little temperature dependence of the masses below T = O.llGeV. While above that 
value, the masses decrease with increasing temperature. The results indicate that the 
hadron-quark phase transition temperature may be T c > O.llGeV for the decuplet 
bayons. 
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I. INTRODUCTION 


The study of hadronic properties under extreme conditions of temperature and density is 
an interesting and important theoretical thesis, being related to present and future heavy- 
ion experiments. Whereas the calculation of hadronic parameters at finite temperature 
and density directly from QCD faces some difficulties. The thermal QCD is a calculable 
theory in large momentum transfer region or at short distance, where the running coupling 
constant is small, so the perturbative method can be used efficiently. However, the coupling 
constant becomes large and the standard perturbation theories fail at the hadronic level. 
Therefore, investigation of hadron properties requires nonperturbative methods. Some 
nonperturbative approaches have been put forward since 1970s, such as lattice QCD, heavy 
quark effective theory (HQET), different quark models, QCD sum rules and so on. In 
addition, it is believed that the hadronic matter undergoes a phase transition to a quark- 
gluon plasma(QGP) at sufficiently high temperature. Whereas the nature of this transition 
is not well understood, so investigation of thermal properties of hadrons is very instructive. 

The QCD sum rule method has been widely used to study the thermal properties of the 
light-light mesons in Refs. jT], 0| , th e heavy-light mesons in Refs. {sUsJ and the heavy-heavy 
mesons in Refs. j6-ll| since it was firstly extended to finite temperature in Ref. 12] by 
Bochkarev and Shaposhnikov. The extension bases on the following assumption: Both the 
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operator product expansion(OPE) technique and notion of quark-hadron duality remain 
valid at finite temperature, but the vacuum condensates must be replace by their thermal 
expectation values. For baryons, Yuji. Koike investigated the octet baryon masses at finite 
temperature, taking into account the contribution of ti+N —> N in the calculation of the 
spectral function 13]. In Ref. [blj], the authors constructed the spectral representation of 
the two-point correlation function of nucleon from the Feynman diagrams. Therefore it is 
interesting to investigate the temperature effect of the baryon masses from the principle 
of QCD. The main aim of this work is to calculate the spectral functions directly from the 
quark propagator at finite temperature. 

The rest of the paper is organized as follow. In section HH we calculate the spectral 
functions to T 8 order, summarize the nonperturbative contributions and obtain the Borel 
transformed sum rules for the masses. Section IIIII is devoted to the numerical analysis and 
a summary is given at the end of this part. 


II. THERMAL QCD SUM RULES OF THE DECUPLET BARYONS 

We begin with the following two-point thermal correlation function: 

IW<Z,T) = i J dx 4 e l<lx (T(J,(x))T u (0)>, (1) 

where J M (x) is the interpolating current for the decuplet baryons. In the above equation, 
the first T is a variable for temperature and the second one stands for time-order product. 
We omit the first one in the rest of the paper in order to avoid confusion. The thermal 
average of any operator O is defined as 

= Tr(exp(-f3H)Q) 

[ ] 7V(exp(-/3if)) ’ 1 ’ 

in which H is the QCD Hamiltonian, — -f stands for the inverse of the temperature, and 
the traces are carried out over the complete set of states. In this paper, we use the Ioffe’s 
currents for the baryons: 

Jp = £abc(u a (x)C'y fl u b (x))u c (x), 

jf = \J^£abc(2(u a (x)C'y ll s b (x))u c (x) + (u a (x)c-f t ,u b (x))s c (x)), 

Jf = \J^£abc(‘2(s a (x)C'y tJ ,u b (x))s c (x) + (s a (x)Cj ll s b (x))u c (x)), 

= £abc(s a (x)C'f fl S b (x))s C (x ), (3) 

where a, b, c are color indices and C is the charge conjugate operator. 

According to the standard procedure of the QCD sum rules, we need to calculate the 
correlator (JT]) in terms of physical particles (hadrons) and in quark-gluon language, and then 
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match the two representations. On the theoretical side, the correlator can be expanded as 
Il(g 2 ) = W ert (q 2 ) + U cond (q 2 ) = [ ds + H cond (q 2 ), (4) 

J s-q 

where p(s) is the spectral density. We calculate the perturbative part of the two-point ther¬ 
mal correlator by using the quark propagator at finite temperature in coordinate spaceflhl]. 

S T (x) = Sq(x) f (nTx), (5) 


where Sq(x) is the quark propagator at T = 0 and / is a temperature modification factor. 
The function f(irTx) satisfies /(0) = 1 at T = 0 and at T ^ 0 

z + l + (z- l)exp(- 2 xj) 


f(z) = zexp(-z) 

which can be expanded at small z as 


[1 - exp(- 2 z )] 5 


/(*) = 1 - 


7z A 

360' 


( 6 ) 


(7) 


For u or d quark, S 0 (x) = ^ 4 , and for s quark, S Q (x) = where m is the 

mass of the s quark. 

After some calculations, it is found that the terms containing T 8 do not contribute to 
the imaginary part of the perturbative side of the two-point thermal correlator. The result 
for the A is 


lulling) = (- 


14 


20(2tt) 3 
1 


15(2tt) 


it T )g^ 4 + 


20(2vr) 


-Q 2 q^, 4 + 


7 


iw' ^ *- mw* T 7 " 7 - 


7 


K T )y v q^ + (- 


-/ + 


7 


-71 


- 4 r 4 ) 7 & 


( 8 ) 


15(2tt) 3 160(2tt) 3 30(2tt) 3 

According to the basic assumption of the extension of the QCD sum rules to finite tem¬ 
perature, we take advantage of the results in Ref. [jllaj] for the condensate part, replacing 
the vacuum condensates by the thermal expectation values, 


Y^cond 


fils' 


(?) = 


■7^(uu)q 2 In {-q 2 ){g^ - ^7/A, + ~\ 


+^{uu) 2 ^{g flu 4- 4 + ^(7 „q v - 7 ^) “ + !»%)} 


< 2^2 H-Q^dsiuaxa^GDig, 


fils' 


Abi,). 


(9) 


Now we turn to the hadronic representation of the thermal two-point correlator ([T]). In 
order to obtain the hadronic representation, a complete set of states A (q,r) are inserted 
into the r.li.s. of Eq. ([I]), 


(r(J„(l))J„(0))=^(T(J„(i))|A(,,r))(A( g ,r)|J„(0)), 


( 10 ) 
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where the summation is made over the spin projection r of the baryon A with momentum 
q, and q 2 = with Ma the mass of A. The coupling of the interpolating current to the 
baryon state is defined as 

(0|J M (0)|A(cj,r)) = A A ^(g), (11) 

in which u^(q) is the wave function of the A baryon, (fa — M^jv^q) = 0, and Aa is a 
constant. Then the hadronic representation of the correlator is obtained as 


^{(J,(i))|A(,,r))(A(,,r)|J40)) 


= — Aa [g^v 4~ g' lv.lv 


+g^M A - -7 /i 7„M A 


.. i / x 2 Wv 4 

+ 3 ( 7 ^ 7^m) 3 M 2 


1(7,- 7„9„)A 


Z A 

2 q^qv 


+ ... 


( 12 ) 


where 


3 M a 

stands for the contributions of continuum states and the interaction between 


the current and particles in the medium 16]. It can be seen that the structures g IMJ and 
g^u 4 contain more information than other structures by comparing the Eqs.Q-Q with 
Eq. (TT2|) . therefore we chose the two structures, g /w and g^ u 4i 1° obtain the sum rule. For 
the structure g^, the spectral function is 


P l(s) — X A M A S(s — M\) + • • • , 


(13) 


and for the structure g lw 4i the spectral representation is 

P 2 (s) = Aa<5(s — M^) + • • • . 


(14) 


According to the quark-hadron duality, the contributions of the continuum states and the 
interaction between current and particles in the medium can be approximated by the OPE 
spectral function. Finally we arrive at 


nf y = I ds 


Pi(s) 


A^M a 


s — q 2 Ml — q A 


+ 


ds 


' S 0 


3^(uu)s 
s — q 2 


W 2 hy = / ds 


P2{s) 


\2 

a a 


+i r 


ds- 


20(2ti-) 


J—_ s 2_14 4rn4 

15(2tt)3 71 1 


'so 


s — q* M\ — q* 7 r 

with so the threshold parameter. On the OPE side, we have 

-1 '~Ujs 

^ 24tt 2 


s — q z 


n 


OPE 


ds 3 tt 2 \ uu ) s + _J_ ln(-g 2 )g s (ua X c^ruG a Xcr ), 


n 


OPE 


'0 


7T 


s — q 


14 yj-4^4 


ds- 


20(2t r) 3 


15(277)' 


s — q z 


+ 




Matching both representations and taking Borel transformation, it is obtained 


(15) 

(16) 


(17) 

(18) 


AUAe-"!!" 3 = 


— ( uu) 
37r 2 


M 4 [ 1 - (1 + -^)e-» /M2 ] + —^e-Ml(uu), (19) 






















5 


x2 - Ml/M 2 _ 1 

A 10(2tt) 4 


M 6 [2 - (2 + ^ + Jj)e^/ M2 ] + ^T 4 A// 2 (1 - e-^/ M2 ) + ^{uu)\ 

( 20 ) 

with the Borel parameter M 2 . In the calculation, we have used the parameterization for 
the mixed condensate 16], 

A Q 


= Mq(uu), 


( 21 ) 


with Mq = 0.8 ± 0.2GeV 2 . 


Equations ffTOl) and (l20l) lead to the mass sum rule for the A baryon 

l^ilabM 4 - ^ M 2 aM 2 
M\ = 3 30 


cM 6 + f tt 4 TW 2 + a 2 ' 


( 22 ) 


where a = -(uu), b = l-{l+jfc)e ~ s °/ M2 , c = 2-(2+j^+^)e~ s °/ M2 , and d = l-e~ s °/ M2 . 

Another way to get the mass sum rule is from Eq. (T20l) and its derivative with respect 
to Aj, which leads to another sum rule for the A mass: 


Af| = ( 


eM 8 + T A bM A 
cM 6 + f 7 i A T A dM 2 + A2zl1^2' 


p/2 


(23) 


with e = 6-(6 + |i + |^ + J 

Similarly, the sum rules for the masses of the £*, 5* and hi baryons can be obtained in 
the same process, taking advantage of the results in Ref. jl7| for the condensates parts. 
The results are, to order m, 


160vr 2 


ML = 


(2a + f)bM A - Ad 2 ( 2 cl + f)M 2 + 5mcM 6 + ^ma 2 + ^ mT A dM 2 


ML = 


J > _ 9 1 j 1 3 ' 3 

cM 6 + fit A T A dM 2 + ^m(4a - f)M 2 + ^{a 2 + 2af) 

2 _ eM 8 + ^T A bM A - ^m(4a - f)M A /2 

S * ~~ cM 6 + f 7 T A T A dM 2 + ^m(4a - f)M 2 + a 2 + 2af) ’ 

( a + 2 f)bM A - 22zdM 0 2 (a + 2 f)M 2 + lOmcM 6 + A2A m (2 a 2 + af ) + m T A dM 2 

cM 6 + f 7 r A T A dM 2 + ^fm{2a + f)M 2 + ^(b 2 + 2af) 


(24) 

(25) 


eM 8 + ^T A bM A - ^fm{2a + f)M A 


:) 


1/2 


m 4 = 


^-2 _ /_ Cihi -r 3 _ 3 _ 

H * cAf 6 + fMTMM 2 + ^m(2a + f)M 2 + A2A (£,2 + 2a f) 
A|A/5m 4 - ^ M 2 fM 2 + 15mcM 6 + 3207r 4 m/ 2 + 56 tt A mT A dM 2 

^71/Tfi I 56 _4nn4 ^7 71 /f2 1 100-9_ £ 1\/T? 1 6407T 4 A9 


_o___ 

cM 6 + f 7 \ A T A dM 2 + 1207r 2 m/Af 2 + A 2 A /2 

, , eM 8 + f tt 4 T 4 6M 4 - 120tt 2 ?u/M 4 x1/ _ 

cA/ 6 + fMT A dM 2 + 120? i 2 mfM 2 + A 2 A /2 

where / = — (ss) and ?n is the mass of the s quark. 


(26) 

(27) 

(28) 
(29) 
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In order to do numerical analysis, we have to know the relations between thermal expec¬ 
tation of operators in our OPE computation and the corresponding vacuum condensate. 
To first order in the pion distribution, the thermal average of an operator O is given by 

3 : 

(0) = (0|0|0)+ Y. f d l kn m l) (^)\0\Ax)). (30) 

,=u.s- / P'O 32 ^ 

Making use of the soft 7r-meson methods, the pion matrix element can be reduced to 
the vacuum expectation value of a double commutator, 

(o) = <oio|o)- T 5^ (si) 

where Q\ — f d 3 xAg(x ) is the axial-vector charge. For light quark case, the result is, 

m = (0\qq\0)(l-^), (32) 

with F n = 0.092GeV being the decay constant of n meson. 


III. NUMERICAL ANALYSIS AND SUMMARY 


Before the numerical analysis of the mass QCD sum rules, we first need to know the 
input parameters of the QCD vacuum condensates. The relation between the condensate 
of u or d quark and the condensate of s quark is (0|ss|0) = 0.8(0|^g|0) 19]. Moreover, we 
use QCD inputs: Mq — 0.8GeV 2 in Ref. 19], (0|tm|0) = — 0.014GeV 3 in Ref. [I2(| and 
m = 0.14GeV in Ref. 21], taking into account the effect of renormalization. 

The sum rules also contain two auxiliary parameters: the Borel parameter M 2 and the 
continuum threshold s 0 - These are not physical quantities, hence the physical observable 
should be approximately insensitive to them. Therefore, we look for working regions of 
these parameters such that the dependence of the mass on these parameters are weak. 
Generally, the continuum threshold s o is related to the square of the first exited state which 
has the same quantum numbers as the concerned hadron, while the Borel parameter M 2 is 
determined by demanding that both the contributions of the higher states and continuum 
are sufficiently suppressed and the contributions coming from higher dimensional operators 
have a good convergence. 

In FigJUwe show the contributions from the excited and continuum states and from the 
higher-dimension operators in the OPE computation for the A case. The lower limit of 
the Borel parameter M 2 is determined by demanding that the contribution of the higher- 
dimension operator is less than 10% of the total contribution and the upper value of M 2 
is the point at which the contribution of the excited and continuum states is 50% of the 
total one. The working window is 2GeV 2 < M 2 < 3.5GeV 2 as we can see from FigJI] In 
addition, the physical quantities should vary weakly with the Borel parameter. Thus the 
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TABLE I: The mass spectra of the decuplet baryons. 


Baryon 

Experimental(Mev) 1201 

This work(GeV) 

A 

1232 ± 3 

1.20 ±0.10 

E* 

1382 ± 2 

1.37 ±0.08 

77 * 

1533 ± 1 

1.50 ±0.09 

n 

2252 ± 9 

2.28 ±0.08 


proper range of M 2 in onr case is chosen as 2.6GeV 2 < M 2 < 3.5GeV 2 . For the parameter 
So, we take values in some range of the square of the first excited state of the baryon A, 
s 0 = 2.4GeV 2 ,s 0 = 2.5GeV 2 ,s 0 = 2.6GeV 2 in our case, and find that its effect can be 
neglected within the current accuracy. Because we don’t know which sum rule is better 
than the other, we analyze the mass sum rules (1221) . (1231) and their average at T = OGeV 
in Fig|2j finding that the average sum rule works well which gives the mass estimation 
M a = 1.20 ± O.lOGeV. 

The same analysis are carried out for the cases of the other three baryons. With the same 
criterion of determining the lower and upper values of the Borel parameter M 2 , the results 
are the following: 2.3GeV 2 < M 2 < 3.5GeV 2 for the E* baryon, 2.5GeV 2 < M 2 < 3.7GeV 2 
for the S* baryon, and 3GeV 2 < M 2 < 4GeV 2 for the baryon. Fig. |3] shows the average 
sum rules for the other three baryons at T = OGeV, from which we get the numerical 
estimates M s * = 1.37±0.08GeV, M=* = 1.50±0.09GeV, and M, n = 2.28±0.08GeV. These 
results show that our sum rules for the masses of the decuplet baryons are meaningful and 
the results agree with the experimental results. In table [J we list the spectra of the decuplet 
baryons and the comparison with the experimental results. 

Fig J4] presents the temperature dependence of the baryons’ masses. It is found that the 
masses nearly vary with temperature below T < O.llGeV, but decrease with increasing 
temperature as T > O.llGeV. Because the hadrons undertake a phase transition from 
hadron states to quark-gluon plasm at sufficiently high temperature, the results show that 
the baryons remain as hadronic states as T < O.llGeV and the critical temperature of 
phase transition to quark-gluon plasma is above T = O.llGeV, ie. T c > O.llGeV. In fact, 
as the temperature is increased, the hadron melts and the width of the hadron increases 
until it becomes infinite at the critical temperature. 

In summary, we study the temperature dependence of masses of the decuplet baryons 
using the thermal QCD sum rules. For the terms containing condensates, we use the results 
obtained by the OPE method and replace the vacuum condensates by the thermal expec¬ 
tation values at finite temperature. The extension of QCD sum rules to finite temperature 
is based on the following assumption: the operator product expansion(OPE) and notion of 
quark-hadron duality remain valid at finite temperature, but the vacuum condensates must 
be replaced by their thermal expectation values. Adopting the quark propagator in the 
coordinate space at finite temperature, we calculate the perturbative part of the thermal 
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(a) (b) 

FIG. 1: The solid line represents the ratio of the contribution of the highest-dimension operator to 
the total contribution; the dashed line stands for the ratio of the contribution of the excited and 
continuum states to the total one. Figs, (a) and (b) represent the g fM/ and the g /JV $ structures, 
respectively. 

two-point correlator to T 1 2 3 4 5 6 7 8 9 order. It is found that the terms containing T 8 order have no 
contribution to the spectral functions and the temperature corrections mainly come from 
T 4 terms. The calculations indicate that the baryon “melts” or the hadron-quark phase 
transition occurs at the temperature T > O.llGeV for the decuplet baryons. 
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(a) (b) 



(c) 


FIG. 3: (a) The average sum rule of X* mass at T = 0 GeV. (b) The average sum rule of H* mass 
at T = OGeV. (c) The average sum rule of Q mass at T = OGeV. 
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FIG. 4: (a) The variations of 
temperature T. (c) The variations of 
with temperature T. 
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